High fidelity one-qubit operations under random telegraph noise 
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We address the problem of implementing high fidelity one-qubit operations subject to time depen- 
dent noise in the qubit energy splitting. We show with explicit numerical results that high fidelity 
bit flips and one-qubit NOT gates may be generated by imposing bounded control fields. For noise 
correlation times shorter than the time for a 7r-pulse, the time optimal 7r-pulse yields the highest 
fidelity. For very long correlation times, fidelity loss is approximately due to systematic error, which 
is efficiently tackled by compensation for off-resonance with a pulse sequence (CORPSE). For in- 
termediate ranges of the noise correlation time we find that short CORPSE, which is less accurate 
than CORPSE in correcting systematic errors, yields higher fidelities. Numerical optimization of 
the pulse sequences using gradient ascent pulse engineering results in noticeable improvement of the 
fidelities for the bit flip and marginal improvement for the NOT gate. 



I. INTRODUCTION 



In physical implementations of quantum computers, 
one of the most challenging tasks is to find an efficient and 
experimentally feasible way to overcome the problems 
caused by undesired interactions between the quantum 
bits, qubits, and their surrounding environment. These 
interactions, which destroy the quantum interference be- 
tween qubit states, lead to errors and loss of fidelity, a 
phenomenon generally referred to as decoherence. 

A variety of methods to fight decoherence have been 
proposed in the literature, including error correcting 
codes P, 0, decoherence free subspace coding 0, 0, 
noiseless subsystem coding dynamical decoupling |y, 
0, quantum feedback control 0, 0, and quantum 
reservoir engineering . Most of these schemes are not 
efficient in making full use of all the physical resources. 
For example, encoding schemes employing decoherence 
free subspaces or error correcting codes store the quan- 
tum information in a specific portion of the whole qubit 
space, or encode several physical qubits into one logical 
qubit. Design and applicability of such codes depends on 
the nature of the decoherence sources and imposes addi- 
tional requirements on encoding and decoding. Dynami- 
cal decoupling schemes possess the attractive feature that 
they require no ancillary qubits, since the interactions 
between the qubits and the environment are effectively 
canceled out by applying external control fields. It has 
also been shown that such decoupling can be realized us- 



ing finite energy soft pulses y|, and that it is possible 
to carry out qubit rotations without disturbing the de- 
coupling process Q. However, dynamical decoupling is 
based on stroboscopic pulsing of the qubit, at a rate sig- 
nificantly faster than the usual characteristic frequency 
of environmental fluctuations. This requires strong con- 
trol fields, which might cause technical problems in the 
laboratory. For example, it has been pointed out that the 
high energy deposition needed for dynamical decoupling 
of nuclear spins is incompatible with the low temperature 
requirement in some qubit implementations |T^|. 

In this paper, we consider the design of fidelity opti- 
mized one-qubit operations in a noisy environment. Mo- 
tivated by experiments in solid state qubits 0, 0] , we 
assume random telegraph noise (RTN) |15| in the qubit 
energy splitting as a phenomenological model for the en- 
vironmental fluctuations. Whereas an ensemble of RTN 
fluctuators models the ubiquitous 1/ / noise in electronic 
circuits [16| , solid state devices on the nanoscale are often 
found to be affected by a single RTN source 0,0] char- 
acterized by its correlation time r c . Although we assume 
that the noise couples only to the qubit energy splitting, 
a similar analysis can be made for errors in the rotation 
angle [HIU. 
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Composite pulses are known to provide an efficient 
way to reduce errors due to systematic off-resonant per- 
turbations, e.g., compensation for off- resonance with a 
pulse sequence (CORPSE) [21|. Here we focus on the 
situation in which the perturbation of a qubit is fluc- 
tuating in time, and seek to suppress the decoherence 
arising from this time-dependent noise by imposition 
of a bounded control field. We are particularly inter- 
ested in the regime where the maximal energy a max pro- 
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vided by the control field is related to the noise corre- 
lation time as T c a max /h ~ 1, since this is an important 
regime for the experimentalists where dynamical decou- 
pling is not applicable. In this regime, we find that a 
novel optimized pulse profile can increase the fidelity of 
quantum operations by up to 30% in comparison to the 
standard composite pulse sequences such as CORPSE 
and short CORPSE (SCORPSE) [UIH. However, we 
provide an analysis over the entire range of correlation 
times, ranging from short noise correlation times satisfy- 
ing T c a max /7i <C 1, where the optimal fidelity is obtained 
by pulsing at minimum time, to the long correlation time 
regime, r c a m&x /h 1, where to a good approximation, 
the fidelity loss is due to systematic time-independent 
errors. The latter case is efficiently tackled by the com- 
posite pulse approach 0, ^ • 

It was stated in Ref. that CORPSE is the short- 
est sequence in the family of composite pulse sequences 
correcting systematic errors as efficiently as possible and 
composing of up to three pulses. Therefore, it was con- 
sidered to be the most useful one. Whereas SCORPSE 
is not as accurate as CORPSE, it is shorter in time. 
Thus it may still be of some interest, depending on the 
physical scenario. We illustrate this fact here by show- 
ing that SCORPSE actually yields higher fidelities than 
CORPSE in the regime of intermediate noise correlation 
time, T c a max /?L ~ 1, for the one-qubit example studied 
here. The optimal performance of CORPSE is in fact 
limited to just the long correlation time regime. We 
also go beyond these composite pulse sequences to ob- 
tain fidelity optimized pulses consisting of large numbers 
of pulse amplitudes that are numerically derived using 
an adaptation of the method of gradient ascent pulse en- 
gineering (GRAPE) H3- Although GRAPE was origi- 
nally developed for finding control pulses in closed quan- 
tum systems, it is utilized here to the determination of 
bounded control pulses in an noisy quantum system. We 
find that numerical optimization of a pulse sequence with 
GRAPE noticeably increases the fidelity of bit flip opera- 
tions, compared with the results of the standard compos- 
ite pulse sequences. In contrast to this improvement for 
a bit flip, only a rather marginal improvement in fidelity 
is found on employing GRAPE to numerically generate a 
fidelity-optimized control pulse sequence for the complete 
one-qubit NOT gate. 

The remainder of this paper is organized as follows. 
In Sec. [n] we characterize the system Hamiltonian, the 
noise model, and the fidelity. Section 1TTT1 introduces the 
pulse sequence generation methods we use for the noisy 
qubit system, and Sec. IIVI presents the results obtained 
for implementation of a state transformation correspond- 
ing to a bit flip and a quantum gate corresponding to the 
one-qubit NOT gate. Finally, Sec. IVl concludes the paper 
with a discussion of extensions and possible generaliza- 
tions. 



II. SYSTEM CHARACTERIZATION 

We consider a single qubit described by the effective 
Hamiltonian 

h= MO + ^WK (i) 

ie{x,y,z} 

where the symbols {a{\ denote the Pauli spin matri- 
ces j24[) {Vi(t)} are the amplitudes of the environmental 
noise, and {cii(i)} are the external control fields. Note 
that the latter are parameterized here by their corre- 
sponding energy amplitudes, rather than the actual phys- 
ical control fields, e.g., electric field amplitudes. We as- 
sume that the strength of the control fields is finite, and 
denote their maximum possible value by a max . To sim- 
plify the discussion, we consider only the case when there 
is no control in the y or z directions, and no noise in the 
x and y directions. Under these assumptions, the Hamil- 
tonian becomes 

H = a(t)a x +r)(t)a z , (2) 

where the control field a :— a x £ [— a max , a max ] and we 
have used the notation rj(t) := rj z (t). 

For RTN, the amplitude of the noise r\ changes ran- 
domly in time between two values — A and A. The quan- 
tity A describes the strength of the noise, and the fre- 
quency of the jumps between —A and A is determined 
by the correlation time r c . Specifically, the probability 
of the noise to jump in an infinitesimal time interval dt 
is given by — . Hence, the probability of no jumps taking 
place in a time interval of length t is 

p (t) = e- t ^. (3) 

In generating sample trajectories of RTN, Eq. J2J can 
also be inverted to yield the sojourn time before a jump 
takes place. Thus we get a sample trajectory of RTN by 
taking random numbers pi € (0,1) and then deriving the 
corresponding jump time instants 

i 

ti=J2~ T oHPj)- (4) 

3=1 

Using the values of these jump times, we can express the 
noise process rj(t) as 

v{t) = (_!)£, eW v (0), (5) 

where 0(f) is the Heaviside step function. 

Since we use an effective Hamiltonian operating solely 
on the qubit rather than treating the full quantum dy- 
namics of both the qubit and the environment, we need 
to average over different noise trajectories in order to ob- 
tain the system dynamics under the influence of RTN. 
Therefore, the dynamics of the system density matrix p 
can be written as 

k=l 
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where po = p(0) is the initial state of the system and 
the operators {U^} refer to unitary time evolution of the 
system under a certain trajectory r]k{t). Formally the 
operator {Uk} is written as 



U k — fe~ i f° dT [ a ( T ) a *+"nk(T)(Jz}/ n ; 



(7) 



where T is the time ordering operator. 

Let p f be the desired final state of the system. Follow- 
ing Ref. [2^ , we define the fidelity function as 



<t>(Pf,Po) = tr{p\p(T)}, 



(8) 



where p(T) is the actual state of the system at the final 
time instant T. Substituting Eq. JBJ into Eq. JSJ, we 
obtain 



N 



<t>{pi,Po) 



J im ij2^{pl u ^oUl}- (9) 

N^oo iv * — ' 
k=l 



Equation (J5J) shows that the fidelity function defined here 
can be viewed as an average over fidelity functions cor- 
responding to individual unitary time developments in 
noiseless quantum systems that are characterized by the 
evolution operators Uk- 

Let us write the initial state of the qubit as 



pa = (I + c x a x + CyCTy + c z a z )/2, 



(10) 



where Cj are real numbers. For implementation of quan- 
tum gates rather than state transformations, we define a 
fidelity function as the average over all pure initial con- 
ditions of the qubit: 



1 

47r ./ c 2 +c 2 +c 2 = 1 



dn^UtpoU^po), 



fin 



where the operator that we desire to implement is de- 
noted by U{ and df2 is an infinitesimal solid angle on the 
Bloch sphere. Simplification of the integral in Eq. 111|) 
yields 



1 1 N 3 

$([/ f ) = _ + li m yy 

{ ' 2 n^oo YIN 

fc=ii=i 



(12) 



III. PULSE SEQUENCES FOR NOISY SYSTEMS 

In this section, we introduce the pulse sequences that 
will be used for suppression of decoherence. We first 
summarize the two composite pulse sequences CORPSE 
and SCORPSE [2l], ^ that were originally designed to 
correct systematic errors in the implementation of one- 
qubit quantum gates. The control fields corresponding 
to the CORPSE pulse sequence are 



a 



ac(t) 



max 1 



for < t' < tt/3 



-<z max , for 7r/3 < t' < 2tt 
flmax, for 2tt < t' < 137r/3, 



(13) 



where t is related to the dimensionless time t' by t' — 
o.ma,xt/h. For the SCORPSE pulse sequence we have the 
control fields 



for 



asc(t) 



-a ma x, for 



< t' < tt/3 
tt/3 < t' < 2tt 
2tt < t' < 7tt/3. 



(14) 



In the absence of noise, the CORPSE and SCORPSE 
pulse sequences generate both the one-qubit NOT gate 
and the bit flip state transformation. 

An alternative to these composite pulse sequences is 
provided by numerical construction of pulse sequences 
optimized for maximum fidelity. Such fidelity-optimized 
sequences may be constructed by an adaptation of the 
GRAPE algorithm which was originally designed to 
steer the dynamics of coupled nuclear spins. No noise 
effects or bounds on control fields are included in the 
original implementation. For full details of the GRAPE 
algorithm for closed quantum systems, see Ref. |23). 

The key feature of the GRAPE algorithm is to ap- 
proximate the continuous pulse shape on a time inter- 
val [0, T] by a function that is constant on n small time 
intervals of length At = T/n, and then to derive the 
corresponding gradients of the fidelity function with re- 
spect to these constant values. Let U™ be the unitary 
time evolution operator corresponding to the time inter- 
val [(m — l)At,mAt] and to the noise trajectory r]k- In 
this interval, the control field is approximated by a con- 
stant, a m . Since the fidelity function <j)(p t , po) is an aver- 
age of the fidelity functions used in Ref. [23j , the gradient 
of <p(p { , po) is obtained as an average of the gradients de- 
rived in Ref. H3. Thus 



Sa m 



iAt 



h n- 



lim 



1 N 



,PTJ}, (15) 



fe=l 



where 



n-l 
k 



TjnTjn- 

PfU k u k 



triUfajUjUkajUl}. and 



Pk 



Tjm+l 
' u k ' 



u k u k ■■■u k p (u k u k ■■■u k y 



(16) 



(17) 



In the case that there exist other control terms {a C k{t)H k } 
in the Hamiltonian, the corresponding gradients can be 
obtained from Eq. (|15fl by substituting a by a k and a x 
by HI 

We note that for A = 0, all the individual RTN trajec- 
tories are identical and consequently the averaging and 
limiting procedures in Eq. (|15(l can be omitted. In this 
case, Eq. l[T5f reduces to the equation for noiseless sys- 
tems derived in Ref. [23j| . 

To derive the gradient of the average fidelity <&, we 
note that 



1 1 3 



(18) 
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Hence, the gradient of Eq. QlSj l can be obtained from 
Eq. JT3J) as 

*fr(l7 f ) = 1 * 5m*jUlaj) 
Sa m 12 Sa m ' ^ ' 

In the GRAPE algorithm, we calculate the gradient of 
the desired fidelity function using Eq. (|T5|) or ifHfl) . and 
update the control fields by moving along the direction of 
the gradient with the restriction a G [— a max , a max ]. This 
procedure results in an optimized pulse sequence for a 
given operation time T. Moreover, the fidelity is also 
optimized with respect to the operation time. 

We note that the pulse sequences yielding the opti- 
mal fidelity for each set of system parameters are not 
unique. In order to find as smooth and as simple se- 
quence as possible, we therefore start from a constant 
control field and use the gradient method to maximize 
the fidelity. To ascertain whether our solution achieves a 
local or the global maximum in fidelity, we repeated the 
procedure for several different, uncorrelated initial val- 
ues of the control field. This resulted in different pulse 
sequences with equal fidelities, suggesting that we have 
indeed found the global maximum, though this cannot be 
conclusively claimed. Thus, when we refer to the results 
of the GRAPE algorithm, we shall describe the corre- 
sponding pulse sequences as optimized rather than opti- 
mal. 



IV. HIGH FIDELITY ONE-QUBIT 
OPERATIONS 

In this section, we present optimized pulse sequences 
implementing high fidelity one-qubit operations that 
were obtained using the GRAPE algorithm, and com- 
pare the results with those from the standard CORPSE 
and SCORPSE pulse sequences. We restrict out atten- 
tion here to two quantum operations on the one-qubit 
system, namely, the state transformation corresponding 
to a bit flip, and the one-qubit NOT gate. 

A. Bit flip 

We consider a one-qubit bit flip, i.e., flipping a one- 
qubit state from one of the two poles of the Bloch sphere 
to the other. This problem may arise, for example, when 
some qubits of a multi-qubit register need to be nipped to 
reach a non-trivial state after a collective initialization. 
The initial and final states can be taken as 

A) = f) and pi = U o) • ( 20 ) 

We first consider the limiting cases: vanishing noise 
correlation time, and infinite noise correlation time. 
Case 1: t c — > 0. In this case, RTN averages out due to 
the well-known phenomenon of motional narrowing pB| 



since the noise changes its sign so rapidly that there is no 
time for the qubit to drift into the direction of the noise 
at any given time. It is therefore not surprising that a 
time optimal 7r-pulsc 

for t G [0,irh/a max ], (21) 

is also fidelity optimal, since it is impossible to correct 
arbitrary fast switching using bounded controls. 
Case 2: r c — * oo. In this limiting case, RTN reduces 
to a constant drift, whereas for large but finite t c the 
drift may be treated as approximately constant. In com- 
parison to a 7r-pulse, pulse sequences such as CORPSE 
and SCORPSE that are specifically designed to correct 
systematic errors will clearly improve the fidelity of the 
desired quantum operation here. 

In Fig. 2] the fidelities obtained from 7r-pulse, 
CORPSE, SCORPSE, and GRAPE pulse sequences for 
a bit flip are plotted as functions of the correlation time 
t c . The noise strength A is chosen to be 0.125 x a max 
in this example. As expected, GRAPE yields the high- 
est fidelities for all values of noise correlation time r c 
since it enforces optimization of the pulse sequence. Note 
that the fidelity curve of GRAPE has a global minimum 
near the correlation time r c « 37i/a max . The existence of 
the minimum is due to the fact that since the GRAPE 
pulse sequences are optimized, they will not only yield 
perfect unit fidelity in the short correlation time limit 
t c — > as a result of motional narrowing, but they will 
also yield unit fidelity in the long correlation time limit 
t c — > oo. The latter argument is true since small system- 
atic errors can be corrected to arbitrary accuracy j2|J 
and provided that the GRAPE algorithm does find the 
global optimal solution. Consequently, there must be a 
minimum at a finite value of r c in the fidelity curve gen- 
erated by GRAPE. The corresponding fidelity curves of 
the CORPSE and the SCORPSE pulse sequences also 
show such minima which result from the fact that these 
sequences are specifically designed to correct the system- 
atic errors. In contrast, the fidelity curve for the 7r-pulse 
is seen to be a monotonically decreasing function of the 
correlation time, reflecting the fact that this pulse can- 
not correct systematic errors. Figure ^ also provides a 
good example of the general result that for intermedi- 
ate noise correlation times SCORPSE is more favorable 
than CORPSE. This is a consequence of the shorter op- 
eration time of SCORPSE, which is more significant at 
finite values of r c than the fact that CORPSE is more 
efficient than SCORPSE in correcting systematic errors. 
This balance between the length of the operation time 
and the accuracy in correcting systematic errors results 
in a cross-over between the fidelity curves of SCORPSE 
and CORPSE at very large correlation times. Hence, 
the CORPSE curve eventually rises above the SCORPSE 
curve for longer t c than shown in Fig. ^ 

As discussed in Sec. IIIII the GRAPE algorithm for 
finding optimal pulse sequences involves optimization of 
the operation time T . However, the RTN used in this 
work does not have any dynamical effect on the initial 
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r c /(7l/o max ) 

FIG. 1: Fidelities (j>(p { ,po) as functions of the noise correla- 
tion time t c for a tv- pulse (Eq. 1)2 If . dotted line), CORPSE 
(Eq. solid line), and SCORPSE (Eq. (0, dash-dotted 

line). The optimized fidelity found using GRAPE is shown as 
the dashed line. The strength of the RTN in this example is 
A = 0.125 



5 10 15 20 25 30 

r c /(7l/o max ) 

FIG. 3: Error e(p f ,p ) = 1 - <t>{Pi,Po) of the GRAPE opti- 
mized pulse sequence, shown as a function of the correlation 
time t c for noise strengths A = 0.25 x <z max (dotted line), 
A = 0.125 x a max (solid line), and A = 0.0625 x a max (dash- 
dotted line). 



density matrix pq. Moreover, the noise is Markovian and 
hence any pulse a'(t) with operation time T' < T may 
be extended to an operation time T without change of 
fidelity by setting 



a(t) 



0, 

a'(t-T 



T 



for < < <T-T 
for T — T' <t <T. 



(22) 



Thus the fidelity is a monotonically increasing function of 
the operation time T. In fact, it is found that the fidelity 
saturates at a maximum value for rather short operation 
times, see for example Fig. |5J 
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FIG. 2: Fidelity (f>(p t , po) of the GRAPE pulse sequence as 
a function of the operation time T, for noise correlation time 
t c = 57i/a max and noise strength A = 0.25 x a max . 

Figure |3| shows the error e(p f ,p ) = 1 — cj>(p { ,p ) of 
the GRAPE optimized result as a function of the cor- 
relation time t c , for several different noise strengths. A 
quadratic dependence of the error on the noise strength, 
e oc A 2 , is observed over the parameter ranges A G 
[a m ax/ 16, flmax/4] and r c G [0,30?L/a m ax]. 



B. One-qubit NOT gate 

In this subsection, we analyze implementations of one- 
qubit NOT gates under RTN. The NOT gate corresponds 
to a 7r-rotation about the x-axis on the Bloch sphere and 
thus also carries out the bit flip considered in Sec. IIV Al 
For quantum gate implementations, we use the fidelity 
function $ defined in Eq. (|11|) . 

The rr-pulse, CORPSE, and SCORPSE sequences are 
specifically designed to implement a NOT gate. How- 
ever, the GRAPE pulse sequences for a bit flip and a 
NOT gate differ since the optimized fidelity functions are 
different for these two operations. Figure 0] shows the fi- 
delities obtained with the 7r-pulse, CORPSE, SCORPSE, 
and GRAPE pulse sequences, as functions of the noise 
correlation time t c . The noise strength is set to be 
A = 0.125 x a max , as in Fig. ^ I n comparison to the 
fidelities for the bit flip shown in Fig. ^ the fidelities for 
the NOT gate derived under the CORPSE and SCORPSE 
pulse sequences are lower, whereas the fidelity under the 
7r-pulsc is higher. Nevertheless, Figs, ^and El show quali- 
tatively the same phenomena, namely, motional narrow- 
ing in the short time correlation limit t c — > and cor- 
rection of systematic errors in the long time correlation 
limit t c — > oo. 

Since we employ the averaged gradient in Eq. 119|) 
which effectively involves three gradients for fixed ini- 
tial conditions, one might conclude that finding the op- 
timized pulse sequences for quantum gates will require 
approximately three times as much computational time 
as for the bit flip. However, as noted above, the GRAPE 
algorithm finds the optimal operation time. This task 
is straightforward in the bit flip case, where as shown 
above, the fidelity is a monotonically increasing function 
of the operation time. For the NOT gate however, the 
optimization is nontrivial. Because of the averaging over 
initial conditions, the optimal fidelity does not necessar- 
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ily increase monotonically with T, as illustrated in Fig. El 
Finding the optimal operation time for a quantum gate 
thus clearly increases the complexity of the problem. 
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FIG. 6: Optimal operation time of pulse sequences obtained 
numerically with the GRAPE algorithm, shown as a function 
of the correlation time r c for noise strength A = 0.125 x a max . 



FIG. 4: Fidelities <6(cr x ) for the NOT gate, shown as functions 
of the correlation time r c for a 7r-pulse (dotted line), CORPSE 
(solid line), and SCORPSE (dash-dotted line). The figure 
also shows the optimized fidelity found using the GRAPE 
algorithm (circles). The RTN strength A is set at 0.125 x 
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FIG. 5: Optimized fidelity $(a x ) obtained from GRAPE, 
shown as a function of the operation time T for noise correla- 
tion time t c = 307i/a m ax and noise strength A = 0.125 x a ma x. 

Additional insight into the efficiency of the GRAPE 
pulse sequence may be obtained by examining the be- 
havior of the optimal operation time as a function of the 
correlation time. As shown in Fig. El the optimal opera- 
tion time of GRAPE increases sharply and approaches 
the time of the SCORPSE pulse sequence at a value 
t c rs I8h/a max . The resulting fidelity also becomes very 
close to that obtained with SCORPSE, see Fig. H ft 
appears from Fig. 0] that errors due to RTN cannot be 
efficiently corrected with bounded controls for correlation 
times shorter than f c , and therefore the optimal opera- 
tion time of GRAPE reduces to that of a 7r-pulse. 



V. CONCLUSIONS 

In this work, we have studied how to perform high fi- 
delity quantum operations on a one-qubit system that is 
subject to random telegraph noise acting on the qubit 
energy splitting. We considered examples of two major 
types of quantum operations, namely, a state transforma- 
tion and a quantum gate. For the state transformation 
we chose a bit flip in which the one-qubit state is flipped 
from the south pole of the Bloch sphere to the north pole, 
whereas for the quantum gate we used the one-qubit NOT 
gate which generates a complete 7r-rotation on the Bloch 
sphere about the x-axis. 

In both cases, we compared the fidelities obtained with 
the standard 7r-pulse, CORPSE, and SCORPSE pulse se- 
quences. The same qualitative phenomena were obtained 
for all three types of pulse sequence. In the limit of van- 
ishing correlation time, motional narrowing occurs and 
implies the 7r-pulse to be the most accurate sequence 
since it is time optimal in implementing the NOT gate 
or a bit flip. On the other hand, the CORPSE sequence 
yields the highest fidelity in the long correlation time 
limit r c — ► oo, since it is designed to efficiently correct 
systematic errors. Over a rather wide intermediate range 
of the correlation time r c , SCORPSE yields the highest 
fidelity among all these three pulse sequences, suggesting 
that it may be a useful approach to suppress environmen- 
tal noise in physical realizations of quantum computers. 

Furthermore, we obtained fidelity optimized pulse se- 
quences using the GRAPE algorithm which was always 
found to yield higher fidelities than the most accurate 
composite pulse sequence. Especially in the bit flip case, 
GRAPE yields noticeably higher fidelities than the com- 
posite pulse sequences. In contrast, GRAPE introduces 
only a rather marginal improvement over the most accu- 
rate composite pulse sequence for the implementation of 
the complete NOT gate. 

The results of this paper provide useful bounds for the 
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implementation of high fidelity one-qubit operations in a 
noisy system without ancillary qubits. Although a simple 
RTN model is used in this paper, we expect that the qual- 
itative dependency of the fidelity on noise strength and 
correlation time will also be present in a general qubit 
system. To investigate the validity of this conjecture, 
one may apply the methods presented here to the study 
of different noise models, e.g., Gaussian noise with a 1// 
spectrum. 

Other extensions of this work include high fidelity con- 
trol of multi-qubit systems. Recent work has addressed 
optimal control of noiseless coupled super conducting 
qubits 27]. For these systems, the environmental noise 
may act on each qubit in either a correlated or uncorre- 
cted fashion, which together with the entanglement of 
the qubits, expands the spectrum of the studies. It is 
an important open question to find a control sequence 
for the inter-qubit coupling term that implements a con- 



trolled NOT gate with high fidelity in the presence of 
noise. To generalize the methods of this paper to noisy 
multi-qubit systems, a reformulation of the equations for 
the fidelity and its gradient is required. 
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